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Context. HD 172555 is a young A7 star belonging to the /? Pictoris Moving Group that harbours a debris disc. The Spitzer IRS spectrum of the 
source showed mid-IR features such as silicates and glassy silica species, indicating the presence of a warm dust component with small grains, 
which places HD 172555 among the small group of debris discs with such properties. The IRS spectrum also shows a possible emission of SiO 
gas. 

Aims. We aim to study the dust distribution in the circumstellar disc of HD 172555 and to asses the presence of gas in the debris disc. 
Metiiods. As part of the GASPS Open Time Key Programme, we obtained Herschel-PACS photometric and spectroscopic observations of the 
source. We analysed PACS observations of HD 172555 and modelled the Spectral Energy Distribution (SED) with a modified blackbody and the 
gas emission with a two-level population model with no coUisional de-excitation. 

Results. We report for the first time the detection of [01] atomic gas emission at 63. 18 yum in the HD 172555 circumstellar disc. We detect excesses 
due to circumstellar dust toward HD 172555 in the three photometric bands of PACS (70, 100, and 160 /jm). We derive a large dust particle mass 
of (4.8 ± 0.6) X lO""* Me and an atomic oxygen mass of 2.5 x lO^^R^ M^, where R in AU is the separation between the star and the inner disc. Thus, 
most of the detected mass of the disc is in the gaseous phase. 
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Circumstellar discs play a key role in understanding how planets 
form and evolve, since it is thought that they originate in gas-rich 
circumstellar discs around young stars, the so-called protoplan- 
etary discs. Around older stars, with ages greater than about 5 to 
10 Myr, another type of disc is seen: dusty, optically thin debris 
discs that are produced by destructive collisions between plan- 
etesimals formed in the protoplanetary disc phase. The younger 
debris discs (~ 5 - 100 Myr) are the likely sites of ongoing ter- 
restrial planet formation, while the older ones (0.1-1 Gyr) cor- 
respond to the "clearing out" phase, when most planetesimals 
left over from planet formation are removed from the system. 

While debris dust has been extensively observed, previous 
studies have provided little unambiguous information on the 
chemical composition of the dust. Most debris discs have fea- 
tureless IR spectra, since the grains are large and cold. However, 
there are a few debris discs that show strong solid-state mid-IR 
spectral features, indicating the presence of small, warm grains. 
These debris discs are probably highly transient and are pro- 
duced during periods of intense coU isional grinding between 
large bodies dLisse et al. 1 12001 l2009l) . Gas can provide addi- 
tional compositional information, but it is rarely detected in de- 
bris discs. Only two - 49 Ceti and HD 21997 - show any trace 
of sub-mm CO emission, indicating that typical debris gas abun- 
dances are low re l ative to those in y ounger protoplanetary discs 
dPent et al.1 120051: iMoor et alj|201 lb . The only debris disc with 
a fairly complete inventory of its gaseous species is /3 Pictoris, 
primarily obtained through UV/optical absorption spectroscopy 



(iRober^e et al.ll2006l) . In general, the debris gas studied so far 
appears to be primarily low density and ionized, because these 
low-density discs are optically thin to dissociating UV photons. 



HD 172555 (iGrav et al.1 l2006l) is one of these rai-e stai's 
with an observed wa rm debris disc. This A7V-type star 
located at 29.2 pc (iGrav et al.1 l2006h belongs to the 
Pictoris Moving Grou p (BPMG), and as such, its age is 
in the range 12 Myr (Zuckerman & SongI l2004l) to 20 Myr 
(Barrado v Navascues et al. 1999). HD 172555 was originally 
found to harbour a bright debris disc through IRAS measure- 
ments of excess thermal emission attributed to circumstellar dust 
grains (Cote 1987). The Spitzer/IRS spectrum revealed strong 
solid-state emission features in the mid-IR, indicating that it con- 
tains a relatively large a mount o f warm dust in a fine //m-sized 
population at T ~ 300 K (iLisse et al. 2009). It has been proposed 
that this disc is the aftermath of a catastrophic collision between 
planetary mass bodies. 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



In this paper, we present Herschel (iPilbratt et al.l 1201 Ol) 
far-IR photometry and spectroscopy of HD 172555, obtained 
as part of the "Gas in Protopolanetary Systems" (GASPS) 
Herschel Open Time K ey Programme with the instrument PACS 
jPoditsch etalj|2010h . In Sec. 2 we describe the PACS spectro- 
scopic and photometric observations. In Sec. 3 we describe the 
main observational and modelling results for the dust and gas 
discs. In Sec. 4, we discuss possible origins for the detected OI 
gas. 



1 



P.Riviere-Marichalar et al.: HD 172555: detection of 63 fim [01] emission in a debris disc 




0.010 



0.001 



1000.0 



Fig. 1. HD 172555 SED. The dashed lin e represents the p hoto- 
spheric contribution fitted using VOSA dBavo et al.ll2008l) with 
Teff = 7800 K, logg = 4.5, Lstai = 7.8 Lq, while the solid one 
represents a photosphere plus modified blackbody model with 
r = 280 K and /3 = 0.2. Blue dots represent Heschel-PACS 
photometric points. 



Table 1. HD 172555 Herschel Space Observatory photometry 
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Statistical 
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2.4 
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3.5 




63.1 63.2 63.3 



63.1 63.2 63.3 



63.1 63.2 63.3 



Fig. 2. [OI] ^Pi -^^ P2 spectral feature. The vertical red lines 
show the positions in the rest frame wavelength for the [OI] 
^'Pi ^''P2 transition at 63.18 ^m.. 



known as spaxels. The observations were centred on the 63.18 
jum [OI] ^Pi ^^P2 transition, within the wavelength range 62.94 
to 64.44 jt/m, with a spectral resolution AjAA of 3570. We also 
observed the star with PACS in pointed chopped range spectro- 
scopic mode on September 10, 2011, with the red channel cen- 
tred around 158 jum, with a spectral resolution of 1265. All spec- 
tra were reduced using HIPE 7 and a modified version of the 
reduction pipeline, including saturated and bad pixel removal, 
chop subtraction, relative spectral response correction, and flat- 
fielding. The spectra were extracted from the central spaxel and 
aperture corrected to account for dispersion in the surrounding 
spaxels. The final, reduced line spectra are shown in Fig.|2] 



2. Observations and data reduction 



3. Results 



2.1. Herschel photometry 

HD 172555 falls into the field of view of our scan maps (obsids: 
1342209059, 1342209060, 1342209061, 1342209062) centred 
on CD-64 1208, an MO-type star with which it forms a wid e 
binary system on a ~ 2000 AU orbit (iFeigelson et al ] l2006h . 
Photometry was acquired with standard minimap configuration, 
with a scan length of 3' and scan angles 70° and 110° for the 
70 and 100 jjm bands. The Herschel-PACS photometer observes 
simultaneously in 70/100 fim and 160yum, so this gave us four 
maps at 160 fim. Photometric observations were reduced using 
HIPE 8. We combined the final images at the same wavelength 
using IRAF (http://iraf.noao.edu/) IMCOMBINE, weighting the 
images with the exposure time. We performed aperture photom- 
etry with an aperture radius of 6" for the 70 and 100 fim bands 
and 12" for the 160 yum band. Photometric errors were computed 
by averaging the standard deviation in nine positions surround- 
ing the source scaled with the square root of the number of pix- 
els inside the aperture. Finally, we applied the aperture correc- 
tions supplied by the PACS team. Photometric fluxes are listed 
in Table[T] The source is point-like in the three PACS bands. The 
PACS photometric points are shown in blue in Fig.[T] 

2.2. Herschel spectroscopy 

The star was observed in chop/nod line spectroscopic mode with 
PACS on March 8, 2011 and September 10, 2011. PACS allows 
imaging of a 47" x 47" field of view, resolved into 5x5 pixels 



3.1. Star properties and evolutionary status 

A compilation of data from the literature was made to model 
the Spectral Energy Distribution (SED, see Table IB. ft . We also 
retrieved IRAC observations from the Spitzer Space Telescope 
archive (AOR 3923456 for 3.4, 4.5 and 5.0 yum bands and 
25202432 for the 8.0 jjm) and extracted aperture photometry 
using MOPEX. The IRAC 3.6 and 4.5 ;um photometry is satu- 
rated. The stellar parameters were computed by fitting Johnson, 
Stromgren, and 2MASS photometry with Kurucz photospheric 
models making use of th e Virtual Observatory SED Analyzer 
(VOSA, Bav o et al.ll2008l) . W e find a value of 7.8 + 0.7 Lq for 
the steUar luminosity (while IWvattet al1l2007i, fin ds 9.5 L(7i) 



Chen et al 



Lisse et al 



and 7800 + 200 K for the effectiv e temperature, whil e 
(2006) finds Teff = 8550 K, and lWvattetaII(l2007h : 
(2009) used Teif = 8000 K. 

Fig. |3] shows the position of the star in the Hertzsprung- 
Russell (HR) diagram. An error of +200 K was adopted for T^s 
and the error bar in the luminosity corresponds to the uncer- 
tainty of 0.2 pc. The evolutionary tracks and isochrones from the 
Yonsei-Yale group ( Yi et al. 2001) were used. Since th e metallic- 
ity of this object is almost solar ([M/H1=0.09 jGrav et al.1 (120061) 1 
the set with Z = 0.02 (solar) was chosen. According to this com- 
parison, we find ^ 1.68 Mq, and logg* - 4.28, which im- 
plies Ri, - 1.56 Rq. Given the position of the star and the un- 
certainties, we can only set a lower limit to the age of ~ 17 Myr 
(Fig. [3]l. The different ages quoted for BPMG members range 
from 12 to 50 Myr, although the latter is an upper limit. 
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Fig. 3. HR diagram for HD 172555. Dashed lines depict the po- 
sition of the theoretical isochrones for different ages (Yi et al. 
1200 ll) while sohd lines depict evolutionary tracks for different 
stellar masses. 

3.2. The dust disc 

As can be seen in Fig. [T] little or no excess is detected at wave- 
lengths shorter than 7 fim, while significant excess is seen at 
longer wavelengths. 

We used a modified blackbody model to fit the IR excess in 
HD 172555: Fcont ^ By x (Ao/A'f where is the Planck func- 
tion, Ao - 13 jjm and /3 -Q for A < Aq. We do n ot use the 870 fim 
LABOCA point reported bv lNilsson et al ](|2009) because the as- 
sociation of the detected 2-cr signal with HD 172555 is ambigu- 
ous. 

The modified blackbody model gives a temperature of 280 + 
9 K, with an opacity index (for A > Aq) /3 = 0.2 + 0.2, while a 
simple blackbody fit gives a higher t emperature of 329 + 6 K, 
but produces a poor fit to the SED. iRebull etaLl (|2008|) fitted 
the HD 172555 IR excess with a blackbody at 310 K, in good 
agreement with our determination. We note that we can fit the 
whole SED with a single modified blackbody, implying that we 
do not see evidence from this data set for a reservoir of colder 
dust producing the far-IR emission. This means that the same 
dust population that produces the mid-IR excess can also pro- 
duce the entire Herschel/PACS excess. 

Using Rin > jR*{jt:f^'^^'^ (lBeckwithetal.lll990 h. we 
computed the minimum radius for the dust distribution. We get 
^in.min = 4.0 + 0.3 AU. The infrared excess computed by inte- 
grating the mo del is Lir/L,. ^ 7.3 x 10"**, slightly higher than the 
value shown in lWvatt et aP (l2007h of 5 x 10 We used the flux 
at 160 yum to estimate the dust mass of the dust disc using Mjust - 
Fy{l60fim)D- /{KyBy(Tdiist)), valid for optically thin discs, where 
Fv(160;um) is the integrated flux density at 160 fim, D is the 
distance to the star (29.2 pc). By can be approximated with the 
Rayleigh-Jeans regime and Ky - 2 x {l.3mm/A'f cm^g"'. The 
final dust mass is (4.8 + 0.6) x 10""* M^. 

3.3. The gas disc 

The [OI] ■'Pi — >"'P2 transition at 63.18 fim was detected in both 
observations, with a signal-to-noise ratio of ~4 in the average 
spectrum. The emission is only detected on the central spaxel, 
indicating that it is centred on the star Line fluxes (Fune) were 
obtained by fitting a Gaussian plus continuum fit to the line and 
are shown in Table |2] In both spectral line observations, the 
line is spectrally unresolved within the uncertainties. The [CII] 
^P3/2 —^^ Pi/2 emission line at 157.74 fim was not detected; in- 



of 63 fim [01] emission in a debris disc 



Table 2. HD 172555 PACS line observations. 



Obs ID 


Center 




Width 




[OI] flux 




Qjm) 




(fim) 




(10 


-"*W/m2) 


1342215649 


63.183 +0.005 


0.025 + 


0.014 


8.2 


±3.0 


1342228417 


63.195 = 


t 0.005 


0.018 ± 


0.005 


9.1 


±3.8 


Average spectrum 


63.188 = 


t 0.004 


0.032 ± 


0.008 


9.2 


±2.4 



stead, we determine a 3-cr upper limit of 2.3 x 10"'^ W/m^ for 
its flux. 

We computed the total [OI] mass using (see AppendixlBb 
Nto, ^ 9.5 X io36 ^..ne(ergs-'cm-^)xDW) , 

|^A„,x(r) 

where N[oi] is the number of oxygen atoms, R is the inner radius, 
Aui is the Einstein A coefficient, and x is the fractional popula- 
tion for the first excited state of the [OI] ^Pi — >^^P2 transition. 
If the emission originated at the minimal possible position of 
flie dust ring at ~ 4.0 AU, the [OI] mass is ~ 0.39 Me. The 
computed mass is huge, so the reservoir for the [OI] gas must 
have been very large. Nevertheless, the total atomic oxygen mass 
strongly depends on the location where the emission comes from 
(as R^). For instance, the [OI] gas mass decreases to between a 
few 10""* Me and a few 10"^ M^ if it is located near the size-, and 
composition-dependent dust sublimation radius (~0. 1-0.5 AU). 
But we note that, even in this case, the [OI] gas mass remains 
large compared to the dust mass. 

4. Origin of tlie gas 

Our HerschelfPACS observations reveal a copious amount of 
atomic oxygen, although the exact mass depends somewhat on 
its unconstrained radial location (Sec. 13. 31 ). This raises the ques- 
tion of its origin. HD 172555 is a rare debris disc. The dust 
is located in the inner system at only a few AU from the star 
(Sec. 13. 21 ). in contrast to typical debris discs composed mainly 
of cold Kuiper belt-like dust at much larger distances from the 
star (e .g. the co-eval BPMG member HD 1 8 1 327. ILebreton et'al] 
I2012h . Here, we speculate about two different origins for the ob- 
served [OI] gas. 

4.1. Recent release during a catastrophic event 

The models of IWvatt et al] (|2007|) argue fliat HD 172555 in- 
frared fractional luminosity is ~ 86 times higher than the maxi- 
mum value expected from steady-state collisional evolution un- 
der their modelling assumptions at an age of 12 Myr. Although a 
different parameter selection could perhaps explain the departure 
from the maximum value, the authors pointed out that a major 
stochastic event may be playing a role. 

Emission featu res in the Sp itzerflRS spectrum of HD 172555 
were modelled by iLisse et al] (12009.) . who suggested the pres- 
ence of glassy silica. They also tentatively attributed faint fea- 
tures to rovibrational transitions of SiO gas. They concluded that 
the dust composition, and the presence of SiO gas (if confirmed), 
point to a violent production process, such as an hyperveloc- 
ity collision between massive planetary bodies. This suggestion 
was additionally suppor ted by the hyperveloc ity collision ex- 
periments conducted by iTakasawa et al.l (1201 ll) . These authors 
found that the ejecta typically had a steeper grain size distri- 
bution than expected for an equilibrium distribution i n a st an- 
dard collisional cascade, as proposed by iLisse et alJ (l2009h to 
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describe the warm dust in HD 172555 (dn/da ~ a "*). This would 
argue for a source of fresh material, which would have had to be 
produced within the las t 0. 1 Myr to be consistent with SiO re- 
condensation. Recently, iPahlevan et alJ (|201 ih showed that sili- 
cate constituents (Si, Fe, Mg and O) can be produced in the form 
of gas during hypervelocity collisions. Therefore, one can spec- 
ulate that the origin of atomic oxygen detected with Herschel 
could be the violent event proposed bv Lisse et al.. C2009i) to ex- 
plain the shape of the Spitzer-\KS spectrum. 

4.2. Gas accumulation over time 

The extremely low sensivity of Ol to the radiation pressure 
from this A-type star (Fernandez et al. 2006) opens the door to 
alternative, perhaps less catastrophic scenarios as for the ori- 
gin of the detected Ol gas. Atomic oxygen could indeed ac- 
cumulate over time at the place where it is produced. For in- 
stance, the considerable collisional activity in the dust ring pre- 
sumably releases a population of poorly bound, or unbound 
grains due to radiation pressure, which can collisionaly vapor- 
ize larger gr ains due to high relative vel ocities. This has been 
proposed by Czechowski & MannI (l2007h to explain the stable 
gas in the co-eval j3 Pictoris system. Assuming that oxygen is 
originally locked into olivine grains, a dust mass of ~ 0.9 
would need to have been vaporized in less than 12 Myr, corre- 
sponding to a mass rate higher than 7 x 10"** M^/yr (10'^ g/s). 
Interestingly, this value com pares to the mass rate derived by 
ICzechowski & Mann (2007) at the peak position of the jS Pictoris 
dust disc (50-100 AU). 

The sublimation of grains close to the star, at ~0. 1-0.5 AU, 
could also contribute to the production of Ol. Althou gh the 
4 AU- radiu s dust ring is collisi on-dominated, models by I WvattI 
(l2005h and iLohne et al] (1201 2|) show that, even in this case, a 
fraction of the dust grains can escape the collisional cascade 
and migrate toward the star by Poynting-Robertson (P-R) drag. 
Silicates extracted from the ring will eventually reach suffi- 
ciently high temperatures as they approach the star and pro- 
duce O, SiO and MgO, and ultimately Si, Mg and O. Assuming 
again that O is locked in olivine grains, a reference dust mass 
of ~ 10"^ Me would need to have been extracted from the 
dust belt by P-R drag to explain the PACS data over an un- 
known time span. A strict minimum supplying mass rate of about 
10""^ Me/yr is obtained assuming a unit efficiency in atomic O 
production by sublimation over the longest possible time span 
(12 Myr), and assuming that the gas has not viscously flown in- 
ward or outward. On the other hand, the P-R drag timescale for 
bound grain s in the 4-AU radius dust ring, fpR > 7500 yr (e.g. 
IWvatil2005l) . yields an upper limit on the present-day P-R drag 
mass loss rate of a few 10"^ M^/yr This implies that this process 
had to last for longer than a few 10"^yr. 

4.3. Discussion 

Different scenarios can qualitatively explain the Ol enrichment 
of the HD 172555 disc. We note that these are not necessar- 
ily mutually exclusive since gas acc umulation could follow the 
giant collision proposed by iLisse et al.l (l2009l) . if it did not oc- 
cur too recently. For instance, it takes about 1 Myr for the de- 
bris created in a ^iant impact to form an axisymmetric disc 
(Ijac kson & Wvatt' 2012), and be consistent with observational 
constraints (Smith etal. 2012). This may leave sufficient time 
to accumulate Ol gas at some place in the disc. Additional 
processes, such as comet evaporation or grain photodesorption 



(IChen et al.ll2007h . also deserve to be mentioned as possible Ol 
providers. In summary, the origin of the detected Ol gas remains 
open, and detailed modelling of the processes discussed here, as 
well as the detection of atomic or ionized Fe, Mg, and Si will 
help to identify the source of Ol gas. 

5. Summary and conclusion 

Our main results can be summarized as follows: 

1 . We report the detection of the [Ol] line in the debris disc of 
HD 172555 with a flux of (9.2 + 2.4) x lO^'^^W/ml This is 
the first unambiguous detection of gas in this system. 

2. We modelled the SED of the object with a Kurucz photo- 
spheric model (T^eff = 7800 K, log^ = 4.5, Lstar = 7.8 Lq) 
plus a modified blackbody to fit the infrared excess at A > 
7 yum. We obtain Tdust - 280 K and /? = 0.2. We derive a dust 
mass of -4.8 X lO^^Me. 

3. Using purely analytical relations and simple assumptions, 
we derive an [Ol] mass of 2.5 x lO^^R^ where R is the 
unconstrained radial location of the gas. 

4. Although the spatial origin of the line is uncertain, we spec- 
ulate about a possible origin in an hypervelocity collision, 
or in gas accumulation over time due to the low radiation 
pressure on atomic oxygen. 

5. HD 172555 provides a valuable window on the processes 
that occur during the early stages of terrestrial planet forma- 
tion, includin g massive collisions like t he one that produced 
Earth's moon (ICanup & Espositolll996l) . 
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Appendix A: HD 172555 photometry compilation 

Table lB.ll shows a compilation of current literature and Herschel photometry for 
HD 172555. 



Appendix B: IHow to derive the oxygen mass 

In the following we explain how to derive the atomic oxygen gas mass from 
prompt emission level population. To estimate the mass of oxygen gas we con- 
sider the excitation of atomic oxygen to its first fine structure level in the absence 
of a collisional partner This situation happens in very low density environments 
such as debris discs. The main mechanism involved is the so-called prompt emis- 
sion and fluorescence. The prompt emission involves the absolution of a photon 
from the star or from the dust at the precise wavelength of the atomic emission, 
and subsequent re-emission. To model the emission, we assume that the ground 
state is the most populated. The population at steady-state for level ^^P2 is given 
by: 



= -ni(Aio + Bio/io) + "oBoi^io- (B.l) 
Assuming that only the first two levels are populated (« = no + "l), 

(B.2) 



Table B.l. Available photometry 



"1 



1 



n 1 +B10/B01 +A10/B01/10 
since goBoi = giBio and Bio = (c"/2/!v')Aio, the fractional population is 



"1 



1 



" 1 + go/gi + (go/gi){2hv^/c2Jio) ' 



(B.3) 



For 01 we have go = 5, gi = 3. Jio is computed using the distance-dilluted stellar 
flux in the Rayleigh-Jeans regime at 63 /im: 



Filter 


^0 


mag 




flux 


C flux 


Ref 


- 


(jum) 


(mag) 


(mag) 


(Jy) 


(Jy) 




Stromgren u 


0.35 


6.354 


0.010 


13.69 


0.13 


6 


Johnson U 


0.37 


5.070 


0.025 


17.63 


0.40 


1 


Stromgren v 


0.41 


5.203 


0.009 


40.49 


0.34 


6 


Tycho B 


0.42 


5.015 


0.014 


38.89 


0.50 


3 


Johnson B 


0.44 


4.790 


0.019 


53.87 


0.94 


1 


Stromgren b 


0.47 


4.891 


0.005 


47.21 


0.22 


6 


Tycho V 


0.53 


4.793 


0.009 


45.51 


0.38 


3 


Johnson V 


0.55 


4.990 


0.015 


38.45 


0.53 


1 


Stromgren y 


0.55 


4.779 


0.003 


45.40 


0.13 


6 


Johnson R 


0.70 


4.887 


0.023 


33.40 


0.71 


1 


Johnson I 


0.90 


4.581 


0.025 


35.74 


0.82 


1 


2MASS J 


1.25 


4.382 


0.26 


28.92 


6.93 


2 


2MASS H 


1.65 


4.251 


0.212 


22.24 


4.34 


2 


2MASS Ks 


2.15 


4.298 


0.031 


12.82 


0.37 


2 


WISE! 


3.35 


4.269 


0.094 


6.069 


0.52 


11 


IRACl 


3.56 


- 


- 


4.949 


0.035 


13 


IRAC2 


4.51 


- 


- 


3.238 


0.038 


13 


WISE2 


4.60 


3.775 


0.059 


5.309 


0.29 


11 


IRAC3 


5.76 


- 


- 


2.214 


0.007 


13 


IRAC4 


7.96 


- 


- 


1.458 


0.004 


13 


AKARI9 


9 


- 


- 


1.451 


0.017 


10 


WISE3 


11.56 


3.615 


0.015 


1.134 


0.016 


11 


Si-5 


11.7 


- 


- 


1.155 


0.002 


8 


WISE4 


22.09 


2.348 


0.020 


0.962 


0.018 


11 


AKARI18 


18 


- 


- 


0.921 


0.020 


10 


Qa 


18.3 






1.094 


0.011 


8 


MIPS24 


23.84 






0.766 


0.030 


5 


MIPS70 


71.42 






0.197 


0.020 


5 


PACS70 


70 






0.191 


0.005 


7 


PACSlOO 


100 






0.089 


0.003 


7 


PACS160 


160 






0.036 


0.002 


7 


LAB OCA 


870 






0.040 


0.01 


9 


SIMBA 


1200 






< 0.026 


0. 


12 



1 (R,Y 2kT, 



Notes. (1): Ijtohnson et alj jl966h : (2): 2 MASS catalogu e 
jSkrutskie et all 120061) : (3): Tvcho-2 catalogue (H0getal. 2000D; 
(4) IRAS faint source catalogue; (5): |^buO et al. (2008); (6) 
Hauck c& Mermilliod (1998); (7): present work; (8): Moerchen et al] 
(2010); (9): Nilsson et al. (2009), see text; (10): .Ishihara et al, (,2010[) : 
(11): WISE Preliminary Release Source Catalog. The WISE magnitude 
limits in the Explanatory Suplement indicate that the 3.4, 4.6 and 12 /im 
photometry of HD 172555 is saturated. However, the WISE pipeline 
performs PSF fitting when dealing with saturated sources and produces 
accurate results for bands 1, 3, and 4. Unfortunately, for band 2 at 4.6 
fim the pipeline does not produce good results for saturated sources 
(K. Stapelfeldt, private communication).; (12): Schiitz et al. (200^; 
(13): photometry reduced by us with MOPEX from IRAC archive data. 
IRACl, IRAC2 and IRAC3 measurements are saturated. 



(B.4) 



where R is the distance to the star and R, the stellar radius. We obtain Jio = 
2.273 X 10"'^^ X (1AU/R(AU))- erg s"' cm"^ fj^-l ^^-1^ fractional pop- 

ulation is thus: X = 8.67 X 10"^ 

The number of oxygen atoms given the flux in cgs is then 

N[oi, = 9.5216 X 103« F^^^^(^rg.-^cm-^)xDHpc^) ^^^^^,^ 

where 9.5216X 10^* combines the conversion from arcsec^ to steradian and from 
AU to parsec. 

The effect of including additional excitation paths is negUgible. We compare 
this with a situation when the emission is produced at 1 AU from the star. For 
example, if we include all main fluorescence pumping paths (-P2 —> '5o — > ^Pi, 
^P2 -> '£)2 -> ^Pi, and ^P2 -> '5o ^ 'D2 -> ^Pi), the final mass of oxygen 
is M([OI]) = 0.0247 Me, compared to M([01]) = 0.025 Me obtained when we 
only include the main excitation path. 
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